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Abstract 

Two of the most fundamental properties of the dark matter particle, the mass and 
the lifetime, are only weakly constrained by the astronomical and cosmological 
evidence of dark matter. We derive in this paper lower limits on the lifetime of 
dark matter particles with masses in the range 10 TeV — 10^^ TeV from the non- 
observation of ultrahigh energy neutrinos in the AMANDA, IceCube, Auger and 
ANITA experiments. For dark matter particles which produce neutrinos in a two 
body or a three body leptonic decay, we find that the dark matter lifetime must be 
longer than 0(10^^ — 10^^) s for masses between 10 TeV and the Grand Unification 
scale. Finally, we also calculate, for concrete particle physics scenarios, the limits 
on the strength of the interactions that induce the dark matter decay. 



1 Introduction 



There is currently mounting evidence for the existence of dark matter (DM) in our 
Universe from various astronomical and cosmological observations [Ij. However, very 
fundamental properties of the dark matter particle such as the mass or the lifetime are 
still fairly unconstrained. An important restriction on dark matter models arises from 
large-scale A^-body simulations [21 [3] which indicate that, in order to reproduce the 
observations of large scale galaxy surveys, the dark matter particle moved very non- 
relativistically at the time of structure formation, namely the dark matter particle is 
"cold" . For dark matter particles which reached thermal equilibrium with the primordial 
plasma, this condition is achieved when the dark matter mass is ^ 1 keV 

provided their relic density today coincides with the dark matter abundance inferred 
from the seven-year WMAP data, fioM = 0.23 [5]. Furthermore, for this class of models, 
the partial wave unitarity of the S matrix imphes the upper limit moM ^ 100 TeV [6] . In 
contrast, for particles which were never in thermal equilibrium, the allowed dark matter 
mass window ranges from 10 /ieV (as for dark matter axions [7]) to the Grand Unification 

scale pismn]. 

Besides, the existing evidence for dark matter does not require the dark matter parti- 
cle to be absolutely stable, as commonly assumed in the literature. Namely, whereas the 
observed longevity of the dark matter particle can be attributed to a symmetry which is 
approximately conserved, this symmetry could be broken by physics at very high ener- 
gies, thus inducing the dark matter decay. This rationale is completely analogous to the 
proton stability, which can be attributed to a baryon number symmetry but which could 
nonetheless be broken by new physics at the grand unification scale, in turn inducing 
the proton decay. Conversely, the study of the dark matter stability provides valuable 
information about the dark matter interactions with the Standard Model particles and 
about the symmetries that ensure the longevity of the dark matter particle. 

In this paper we will derive limits on the dark matter lifetime focusing on the case 
where the dark matter particle is superheavy, complementing recent analyses studying 
the stability of dark matter particles with masses in the range 1 GeV-10 TeV [HI [111 
fT3l [HI Uni Uni HZl HH HSl EUJ . The decay of superheavy dark matter particles has been 
considered in the past as a source of high energy cosmic rays [2T | [22 | [23 t [2H [25l[26l[27l[28] . 
In the so-called "top-down" models, the decay of the dark matter particles into partons 
produces a flux of cosmic protons, which is accompanied by a flux of gamma-rays and 
neutrinos. Neutrinos have the property that they travel in the Universe without suffering 
an appreciable attenuation, making them a very suitable messenger to constrain the 
lifetime of superheavy dark matter particles; the possibility of detecting of neutrinos in 
"top-down" models has been discussed in [201 EZl ESI EI] • In this paper we will concentrate 
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on scenarios where the dark matter particle decays leptonically, producing a flux of 
high energy neutrinos. To derive bounds on the hfetime we will use the limits recently 
published by various experiments on the high energy neutrino flux and which cover a wide 
range of energies: AMANDA between 16TeV and 2.5 x lO'^TeV [32], IceCube between 
340 TeV and 6.3 x 10^ TeV |33l[3l], Auger between 10^ TeV and 10^ TeV p5] and ANITA 
between 10^ TeV and 3.2 x 10^^ TeV [M]- All these experiments probe different energy 
windows and perfectly complement each other, allowing us to derive stringent limits on 
the lifetime of dark matter particles with masses ranging between 10 TeV and the Grand 
Uniflcation scale. 

The paper is organized as follows: in sec. [2] we discuss the expected neutrino flux 
from both two body and three body decays of dark matter particles, taking into account 
the opacity of the Universe to ultra high-energy neutrinos. Using the expected neutrino 
flux, we derive in sec. |3] lower limits on the dark matter lifetime from the data obtained 
by the AMANDA, IceCube, Auger and ANITA experiments. In sec. |4] we translate the 
derived limits on the dark matter lifetime into limits on the coupling constants of effective 
Lagrangians describing concrete scenarios of dark matter decay. Lastly, in sec. |5] we 
present our conclusions. 



2 Neutrino flux from dark matter decay 

The neutrino flux from dark matter decay receives two contributions. The flrst one stems 
from the decay of dark matter particles in the Milky Way halo and leads to a differential 
flux given by: 

oo 

"^^'^'^ b) = ^ / ds phaio[r(., /, b)] , (2.1) 



dE^ Ati moM tdm dE, 





where phaio('") is the density proflle of dark matter particles in our Galaxy as a function 
of the distance from the Galactic center, r, and dNy/dEy is the energy spectrum of 
neutrinos produced in the decay of a dark matter particle. The neutrino flux received at 
Earth depends on the Galactic coordinates, longitude / and latitude 6, and is given by a 
line-of-sight integral over the parameter s, which is related to r by 



r(s, /, h) = sJs'^ + R%- 2sRq cos b cos I , (2.2) 



where Rq = 8.5 kpc is the distance of the Sun to the Galactic center. 

The average flux over the full sky can then be straightforwardly calculated, the result 
being 

C^-'halo p. dN^ 

-dE7 ~ ^''^''Jk ■ ^^-^^ 



For our numerical analysis we will adopt a Navarro-Frenk- White density profile 

Phaio(r) , , , (2.4) 

r/rc(l + r /rcY 

where r is the distance to the Galactic center, ~ 20kpc is the critical radius 
and ph — 0.33 GeV cm~^, which yields a dark matter density at the Solar System 
Pq = 0.39GeVcm~^ [38]. For this choice of parameters, 

Dhaio = 1.7 X 10-« (—] (—] cm-2 s-i sr-^ . (2.5) 

In addition to the neutrino flux from the decay of dark matter particles in the 
Milky Way halo, there is a second contribution stemming from the decay of dark matter 
particles at cosmological distances, which produces a perfectly isotropic diffuse neutrino 
flux. The differential flux with respect to the received neutrino energy is given by: 

oo 

dJeg _ ^dmPc f , 1 dNi 



dE^ 47rmDMTDM 







where H{z) = /Joa/^a + ^m(l + z)-^ is the Hubble expansion rate as a function of 
redshift z and pc = 5.5 x lO^^GeV/ cm^ denotes the critical density of the Universe. 
Throughout this work we assume a ACDM cosmology with parameters Q\ = 0.73, 
= 0.27, ^DM = 0.23 and h = i/o/lOOkm s'^ Mpc"^ = 0.70, as derived from the 
seven-year WMAP data [5]. Lastly, Su{Ei,,z) is the neutrino opacity of the Universe, 
which was calculated in [291 EH], assuming that the neutrinos are massless, for an age 
of the Universe to = 6.5 Gyr. With the most recent determination of the age of the 
Universe, to = 13.76 Gyr [3], the neutrino opacity reads: 

7.4 X 10-1^(1 + z)V2(^^/TeV), for 1 « z < ,^ ^, 

s^{E^,z) = < (2.7) 



1.7 X 10-1^(1 + zfiEj TeV), for z > 2. 



eq 



where Zeq ~ 10^ is the redshift of the equality of matter and radiation in the Universe. 

In this paper we will analyze the following scenarios. First, we will study the limits 
on the lifetime of a scalar dark matter particle which decays into a neutrino and an 
antineutrino, 0dm — ^ ^i^- In this case, the energy spectrum reads: 



Our second scenario consists of a Majorana dark matter particle which decays into a 
photon and a (anti)neutrino, ip — )■ 72^(2?). Due to the Majorana nature of the dark matter 
particle, the branching ratio for each decay channel is, at tree level, equal to 50%. In this 
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Figure 1: Contributions to the neutrino flux at the Earth from two body and three body 
decay of a dark matter particle with m^M = 10 TeV and r^M = 10^® s. 



scenario the energy spectrum of the (anti)neutrinos is also given by Eq. (2.8) and displays 
a very sharp line feature. To evaluate the limits on models which produce neutrinos with 
a softer spectrum we will also study a scenario consisting of dark matter particles which 
decay into an electron-positron pair and a (anti)neutrino if) — )■ e+e^z/(z/). As before, the 
branching fraction for each decay channel is ~ 50%. In this case the energy spectrum 
of the neutrinos is fairly model dependent, namely it depends on whether the decay is 
mediated by a heavy charged scalar or a heavy charged vector and on the chirality of 
the Standard Model fermions to which the dark matter particle couples [12]. When the 
dark matter particle couples just to the left-handed electrons we find 

dN^ 2 / 2E^y ^ 



(2.9) 



dEi, \ mDM / 

both in the case when the decay is mediated by a scalar or by a vector. 

We show in Fig. [T] the various contributions to the differential neutrino flux for the 
two body and for the three body dark matter decays, Eqs. ( 2.8[2.9 ) respectively. In this 
figure we assumed ttidm = 10 TeV and tdm = 10^^ s. As can be seen from Fig. [l| the 
maximum energy of the neutrino is mDM/2 and the vertical black solid line represents 
the sharp line in the neutrino spectrum coming from the two body decay — )■ z/z/ of 
dark matter particles in the galactic halo. The red solid line shows the extragalactic 
contribution for the two body dark matter decay (see Eq. (2.6)). The dashed (blue) and 
dotted (green) curves show respectively the galactic and extragalactic contributions for 
the three body decay if) — )■ e^e~v. It is apparent from the plot that the neutrino flux 
at the energies near to mDM/2 is dominated by the halo contributions while at lower 



energies by the extragalactic contribution. This feature exists for all the dark matter 
masses and, as we will see in the next section, will play a role in extracting the lower 
limit on the lifetime of superheavy dark matter particles. 

The initial flavor composition of the neutrino flux from decaying dark matter at the 
production point 7° : J° : J° depends on the details of the dark matter model. After 
production, neutrinos travel long distances to the Earth (of the order ~ kpc for the 
Galactic decaying dark matter particles and > Mpc for the extragalactic ones). Due 
to these long distances, the neutrinos arrive at the Earth decoherently, such that the 
flavor composition of the neutrino flux at the Earth is (J®, J®, J®)'^ = P {J^, J^, Jr)^, 
where the elements of the symmetric oscillation probability matrix P are given by Pap = 
J2i \ Uai\^\Ui3i\'^- Assuming the best-flt values for the mixing parameters |10] including the 
updated result on ^^13 from the recent measurements (sin^ 26'i3 
matrix P is given by 



0.1 PT]), the symmetric 
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(2.10) 



Comparing the values of the elements of matrix P, it follows that even if the dark 
matter decay produces neutrinos with just one flavor, similar fluxes are received at the 
Earth for all the three flavors. Therefore, in our numerical results we will assume for 
simplicity that the flavor composition of the neutrinos are in the ratios J® : J® : J® = 
1:1:1. Any deviation from this assumption leads to at most a factor of two difference 
in our limits. 



3 Limits on the dark matter lifetime 



In this section we will present the limits on the dark matter lifetime from various ex- 
periments including AMANDA [32], IceCube-22 [33], IceCube-40 [31], Auger [35j and 
ANITA [SH]- However before discussing each experiment, let us flrst describe the method 
that was employed to derive the limits. Generally, for each of these experiments, it is 
possible to deflne an "effective area" A"^ for the detection of Ua + Pa {oi = e,^,T). By 
deflnition, the effective area is the hypothetical area of the experiment with 100% effi- 
ciency which can be obtained after developing appropriate cuts for the background-signal 
discrimination. The expected number of events is given by 



exp 



dX 



A:^{E,)dE, 



(3.1) 
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where dJ®^_^_p^/dEj, is the energy spectrum of z/q + z/q flux at the surface of Earth, T is the 
hve time of the experiment, AQ is the sohd angle acceptance of the experiment, and 
and E'^^^ are, respectively, the minimum and maximum of the neutrino energy considered 
in the data analysis. The upper limit on the number of neutrinos from the diffuse flux 
depends on the number of signal candidates A^sig (after performing appropriate cuts) and 
the expected number of background events A^bg- For a given A^sig and A^bg, the Bayesian 
upper limit A^umit on ^cxp at q% confldence level can be extracted from the following 
equation 

where L(A^sigl^) is the likelihood function for Poisson distributed variable A^sig which is 
given by 

L(Ar.,|Ar)= ^^ + y% -("^"^.). (3.3) 

-'^sig- 



After extracting the value of A^umit from Eq. (3.2), the upper limit on the diffuse flux of 
neutrinos at a certain confldence level can be derived from A^exp < A'^iimit . 

Various experiments have recently searched for a diffuse flux of ultra-high energy 
neutrinos. The AMANDA experiment has searched for a diffuse flux of muon neutrinos 
in the energy range 16 TeV < By < 2.5 x 10^ TeV, flnding A'gig = 7 signal events while the 
expected number of background events is A^bg = 6, hence giving a 90% CL. limit on the 
total number of neutrinos in that energy window A^umit = 5.4 [32]. Besides, the IceCube 
collaboration has undertaken a dedicated analysis searching for ultra-high energy tau 
neutrinos with the 22-strings conflguration, without discriminating between electron and 
muon neutrinos The energy range of the analysis is 340 TeV < < 2 x 10^ TeV 
and the expected number of background events is A^bg = 0.6. In this analysis the number 
of candidates of signal events is A^sig = 3, which results in A^umit = 6.1 at 90% CL. An 
extended analysis using the 40-strings conflguration yielded A'sig = in the energy range 
2 X 10^ TeV < Ey < 6.3 x 10^ TeV while A^bg = 0.1 [MJ, which translates into A/'umit = 
2.3. Furthermore, the Auger experiment has covered an energy window with extends to 
higher energies, 10^ TeV < E^ < 10« TeV. In this case, A^'bg = and A^'sig = [35], 
which results in A^umit = 2.3. Lastly, in the extremely high energy region, the ANITA 
experiment has covered the window 10® TeV < E^, < 3.2 x 10^^ TeV, flnding A"sig = 1 
while A^bg = 0.97 [3B], giving A^umit = 3.3. Table [l] summarizes the energy window covered 
by each experiment and the relevant experimental results for our analysis. Also, using 
the information in Table [T| it is possible to obtain the upper limit on the astrophysical 
diffuse flux of neutrinos. Assuming the spectrum kE'"^ for the diffuse flux, we obtained 
upper limits (7.4 x 10-^ 1.6 x 10"^ 3.6 x 10-^ 1.7 x 10"^ 1.3 x 10"^) GeV cm^^ g-^ sr'^ 
on the normalization factor k, respectively for the experiments AMANDA, IC-22, IC-40, 
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Auger and ANITA; which matched the reported hmits by each experiment. 

Table 1: Summary of the experiments used in deriving the hmit on dark matter hfetime, 
with their corresponding energy window of sensitivity, the number of background events, 
A^bg, the number of signal events, Ngig, and the upper limit on neutrino events at 90% 

C.L., A^limit- 







^bg 


' Slg 


A^limit 


AMANDA 


16 - 2.5 X 10^ 


6 


7 


5.4 


IceCube-22 


340 - 2 X 10^ 


0.6 


3 


6.1 


IceCube-40 


2 X 10^ - 6.3 X 10*^ 


0.1 





2.3 


Auger 


10^ - 10*^ 








2.3 


ANITA 


10^ - 3.2 X 10" 


0.97 


1 


3.3 



For each of these experiments, we derive 90% CL. limits on the dark matter lifetime 
as a function of the dark matter mass requiring that the number of expected events 
from dark matter decay does not exceed Nn^n, cf. Table [T} To this end, we use the 
"effective area" reported by each of these experiments and calculate the total number 
of events from Eq.(3.1), taking the neutrino fluxes in Eqs. ( 2.3[2.6 ) for the two and the 
three body decay spectra, Eqs. ( 2.8[2.9 ). Fig. |2] shows the 90% C. L. limit on tdm for 
the scenario where the dark matter particle is a scalar which decays into a neutrino 
and an antineutrino, — )■ z/z/. The horizontal lines at 4.3 x 10^^ s and 2.2 x 10^^ s 
show respectively the age of Universe and the limit derived from the position of the first 
peak of the cosmic microwave background [12] • The excluded regions from left to right 
correspond, respectively, to the experiments Super-Kamiokande, IceCube-22 low energy, 
AMANDA, IceCube-22, IceCube-40, Auger and ANITA. The Super-Kamiokande limit 
was taken from [16], while the "IC-22 low energy" excluded gray region from [l3], which 
performed a dedicated search for the two body decay of dark matter particles in the 
galactic halo employing the IC-22 configuration; all other limits are new, to the best 
of our knowledge. The limits derived for each experiment consist of two parts: 1) a 
peak-shaped region which stems from the galactic halo contribution to the neutrino flux 
and which corresponds to dark matter masses in the range ii^™" < m£)M/2 < E^^^, 
namely to dark matter scenarios where the produced neutrino has an energy falling in 
the window to which the experiment is sensitive; 2) a high-energy tail, which corresponds 
to dark matter masses mDM/2 > E^'^'^, such that only the neutrinos produced in decays 
at large redshifts have energies in the window of sensitivity of the experiment. For each 
experiment in Fig. |2| the solid line in the tail part shows the limit neglecting the opacity 
of the Universe, Sjy, while the shaded area shows the effect of the opacity. As apparent 
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Figure 2: Limits (90% C. L.) on the dark matter lifetime tdm from various experiments, 
assuming two body decay (f) — ?■ vv. The horizontal lines at 4.3 x 10^^ s and 2.2 x 10^^ s 
show respectively the age of Universe and the limit derived in from the position 
of the first peak of the cosmic microwave background. The excluded regions from left 
to right respectively correspond to the experiments: Super-Kamiokande taken from |16j . 
IceCube-22 low energy taken from |13], AMANDA, IceCube-22, IceCube-40, Auger and 
ANITA. 

from the plot, the absorption of neutrinos only plays a role in the extremely high energy 
range and effectively restricts the capabilities of experiments to limit the lifetime of dark 
matter particles with mass around the Planck scale. 

On the other hand. Fig. [3] shows the 90% C. L. limit on the dark matter lifetime 
Tdm assuming the three body decay ip — )■ e'^e^v{v). It is noticeable that, due to the 
larger extragalactic flux at low energies for heavy dark matter particles, the limits from 
the tail part in Fig. [3] are stronger than the corresponding limits in Fig. [2] Furthermore, 
for a given experiment, the limits in Fig. |3]for dark matter masses in the range < 
^dm/2 < E^^^ are approximately a factor of two smaller than those in Fig. [21 This 
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Figure 3: Same as Figure [2| assuming three body decay ijj — )■ e+e v. 



can be easily understood calculating analytically, from Eq. (3.1), the lower limits on 
the two body decay and three body decay lifetimes, t^^^"'^^ and Ti;~|J°'^^. Taking into 
account that the neutrinos from the two-body decay of dark matter particles in the halo 
are monoenergetic and neglecting the contribution from the cosmological decays, it is 
straightforward to calculate t^^^""^^ as the values of tdm which saturate the following 
inequality 

iVexp - ^ Y: ^e-ff (^) < iVlim. , (3-4) 

"^dmTdm 3 ^ V 2 / 



where d is the numerical factor in Eq. (2.5), which is equal to 1.7 x 10~^ cm~^ s~^ sr~^ 
for rriDu = 1 TeV and tdm = 10^^ s, and the factor 1/3 comes from the assumption 
of democratic flavor composition of neutrinos at the Earth J® : J® : J® = 1 : 1 : 1. 
Conversely, using Eq. (3.4) to cast the effective area in terms of Tii~|j°'^'^(^DM), it is easy 
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to show that the hmit on the three body decay hfetime t^^^"'^^ can be written as 

^max 

rLr^(mDM) J ^ (2^^) W '^(2i?.) dE^ , (3.5) 



mm 



where dN,y/ dE^ is given in Eq. ( 2.9 ) and the factor 1/2 in front of the integral comes from 



the fact that the number of neutrinos plus antineutrinos per three body decay is half of 



the number of neutrinos per two body decay. From Eq. (3.5) it can be checked that, for 
a given experiment, when ttidm — 2E^'^'^ the relation Tjj~|J°'^-^ ~ ''"iimit°'^^/2 approximately 
holds, as almost all the spectrum of neutrinos lies inside the integration range of the 
experiment. However, for < 1E^^^ it follows that t^^^°^^ ^ ''"limit"'^'*'/^ since some 
part of spectrum lies outside the integration range of the experiment. This procedure 
has been employed in Fig. |3] to derive the "IC-22 low energy" limits from the published 
lower bounds on the lifetime for the decay — )■ i/z/ [13] ; for all other limits, the effective 
areas provided by the experiments were used. 



4 Implications for decaying dark matter models 

In this section we will translate the upper limits on the dark matter lifetime derived in 
the previous section into limits on the couplings of the effective Lagrangian which induce 
the dark matter decay. Let us first discuss the scenario where the dark matter particle is 
a scalar. In this case the effective Lagrangian that induces the decay 0dm ~^ contains 
a dimension 4 operator: 

- £cff = 1/0dmz^z^ + h.c. (4.1) 

where y is the couphng constant. The decay rate can be straightforwardly calculated, 
the result being: 

|y|2 

r(0DM vv) = — moM . (4.2) 
Stt 

We show in Fig. |4] the upper limits on the coupling y as a function of the dark matter 
mass from the non-observation of neutrino events in the experiments discussed in sec. [3j 
The solid line shows the upper limit obtained in this paper, while the dotted line shows 
the upper limit coming from requiring a dark matter lifetime longer than the age of 
Universe and the dashed line shows the limit from the CMB data |12]. 

On the other hand, when the dark matter is a Majorana fermion, the effective La- 
grangian that induces the decay ip ^ reads [S] 

C = ^^DM (Tapifi + e75) z/F"^ + h.c. (4.3) 



11 




wdm (TeV) 



Figure 4: Upper limit (90% C. L.) on the coupling y which induces the two body decay 
— )■ z/z/ (see Eq. (4.2)). The dotted and dashed lines show respectively the limit from 
age of Universe and CMB. The solid line shows the limit obtained in this paper. 



where F""^ is the electromagnetic field strength tensor, while /i and e are the magnetic 
and electric transition moments, respectively. When the dark matter and the neutrino 
have the same CP parities, the magnetic transition moment vanishes, while when they 
have oppositte CP parities, the electric transition moment vanishes. In either case, the 
decay rate can be cast as: 

r(z^DM ^ vi) = ^-^mU , (4.4) 



where we have defined an effective neutrino magnetic moment, \jj,es\ = \/\JA 



In specific models, the effective Lagrangian Eq.(4.3) is generated by loop effects. 
Concretely, when the decay is mediated by a heavy scalar circulating in the loop, the 
effective neutrino magnetic moment can be conveniently parametrized as [T9] : 



l^^^l = 64vr2M| ' ^^-^^ 

where My, is the mass of the scalar particle in the loop and ^efi is a combination of the 
couplings of the dark matter particle and the neutrino to the heavy scalar in the loop. 
Assuming = Mp, where Mp is the Planck mass, we show in Fig. |5] upper limits on 
the coupling ^^cfr! when < Mp, the limits scale as (M^/Mp)^. 

Lastly, the width of the three body decay of a Majorana dark matter particle into 
an electron-positron pair and a neutrino ipBM e^e~v is fairly model dependent, as 
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Figure 5: Upper limit (90% C. L.) on the parameter ^eff which induces the radiative two 
body decay ip ^ wj (see Eqs. ( 4.4[4.5 )), assuming Ms = Mp. The dotted and dashed 
hues show respectively the hmit from the age of Universe and CMB. The solid line shows 
the limit obtained in this paper. 



discussed in [TSj . Let us first assume that the dark matter decay is mediated by a heavy 
charged scalar E. The Lagrangian reads: 

C^s = -V-DM [Xe^PL + Xe^Pn] c - uXf^PRC + h.c. , (4.6) 

where Pl = (1 — 7^)/2 and Pr = (1 + 7^)/2 are the left- and right-handed chirality 
projectors, respectively. Then, the dark matter decay width is given by: 

I \ |4 ™5 

r(^OM^e-e -) = ^^^ (4.7) 

where |Aeff|'^ = (|A^p + |Ag^p)|A^p. Fig. 6] shows the limit on Agfr as function of ttidm 
assuming Ms = Mp; when M^ < Mp the limits scale as (M^/Mp)'^. 

In contrast, if the decay is mediated by a charged vector, the effective Lagrangian 
read^ 

C = -^BuY [X^^Pl + Af^Pp] e - u^Xf^Pne + h.c. . (4.8) 
In this case the decay width is: 

r„„„^e-e-.)^M;!|l (4,) 



^We assume here that the decay is dominated by the charged-current interaction; in more generality 
the decay could also be mediated by a neutral current interaction. 
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Figure 6: Upper limit (90% C. L.) on the parameter Acfr which induces the three body 
decay ip — )■ e'^e~u (see Eq. (4.7)), assuming Ms = Mp. The dotted and dashed hues 
show respectively the limit from age of Universe and CMB. The solid line shows the limit 
obtained in this paper. 



where, again, |A 



eff 



2 -L I I 



)|A^p. Note that the decay width for the decay 
mediated by a vector is a factor of four larger than the one mediated by a scalar, assuming 
Mv = Ms in Eqs. and K9\). 



5 Conclusions 



We have derived lower bounds on the lifetime of dark matter particles with masses in the 
range 10 TeV < rriDM < 10^^ TeV from the non-observation of high energy neutrinos in 
the experiments AMANDA, IceCube, Auger and ANITA. We have analyzed two scenarios 
where the dark matter decay produces monoenergetic neutrinos, namely a scalar dark 
matter particle which decays into a neutrino and an antineutrino and a fermionic dark 
matter particle which decays into a neutrino and a photon, as well as a scenario where 
the dark matter decay produces a neutrino flux with a softer spectrum, concretely a 
fermionic dark matter particle which decays into a neutrino and an electron-positron 
pair. We have found that, for dark matter masses between ~ 10 TeV and the Grand 
Unification scale, the lifetime is constrained to be larger than O{10^^ — 10^®) s, which 
is nine to eleven orders of magnitude longer than the age of the Universe, and seven to 
nine orders of magnitude stronger than the limit from the cosmic microwave background. 
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We have also translated the limits on the lifetime into limits on the parameters of the 
Lagrangian of our three scenarios. Especially for large masses, the limits on the couplings 
are very strong and hint to the necessity of a symmetry ensuring the stability of the dark 
matter particle. 
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